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The bell prover is widely used for gas 
flowrate measurements by timing a 
known stroke of the bell as it rises, pre- 
sumably with a constant speed, from a 
bath of sealing liquid. A differential 
equation for the bell motion (from the 
previous paper of the same title) is used 
together with the previous equations of 
motion for the gas and sealing liquid 
modified here to provide the basis for a 
computer simulation of the prover per- 
formance. Examples of the results of the 
computations show substantial fluctua- 
tions in all of the motions and modified 
measurement procedures for improved 



accuracy are discussed. Proposed modi- 
fications to the prover itself are shown 
by the computations to dampen the fluc- 
tuations and to improve the measure- 
ment accuracy. Only a limited number 
of changes of the prover design and of 
the initial conditions were researched 
for their effect. 
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1. Introduction 



The bell prover of figure 1(a) is essentially an 
ordinary prover which has been modified with 
parts added in an attempt to provide the basis for 
an improvement in gasflow measurement accuracy. 
A porous plug, a heat exchanger, loose metallic- 
wool (throughout the pipe volume Fb), and thermal 
insulation are installed in the inlet pipe. These are 
suggested to help attenuate the fluctuations of pres- 
sure that are generated in the bell and also to assure 
equality of temperature of the bell and of the gas 
input. A second essential addition consists of two 
annular orifice plates mounted in the two sealing 
liquid spaces between the bell and its tank. It's ex- 
pected that the constricted liquid flow area there 
will help to attenuate the liquid motion which is 
shown in the subsequent computations herein to be 
a major source of flow measurement inaccuracy. 



Three other additions of lesser importance are 
shown in figure 1(a) and these are suggested for use 
with both the modified and ordinary prover sys- 
tems. The first is the utilization of a sonic flow noz- 
zle to be placed somewhere in the bottom of the 
inlet pipe. It will serve to prevent fluctuations in 
the gas motion that are generated in the bell from 
feeding back to the meter thereby causing uncer- 
tainties in its calibration. A second benefit is that it 
eliminates the need to include, perhaps, a different 
pipe volume with the provers tare volume Vi each 
time a different meter is installed for calibration. A 
third benefit arises in that the computations that are 
to follow here will be more hkely to represent ac- 
tual events for this restricted system as compared 
to the usual large complicated piping system. 
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Figure 1. (a) Schematic drawing of the bell prover. (b) Details of liquid seal and its velocity 
profile. 



More than one nozzle (size) may be needed to be 
able to cover the full range of flow of any one 
prover. Other basic requirements are that, at the 
flow rate of choice, the gas velocity in the nozzle 
throat be equal to the velocity of sound there, and 
that the gas at the entry point to bell be at the 
temperature of the bell. Another concept may 
prove useful here. Simultaneous measurements of 
temperature and pressure can be made at the sonic 
nozzles to calibrate them along with the meter. 
These calibrated nozzles then should prove useful 
as secondary flow standards for calibration of me- 
ters without the necessity of always exercising the 
bell prover. A second addition to the ordinary 
prover should include a thermocouple placed at 
the exit of the inlet pipe so that it can provide a 
confirmation that the gas and bell are at equal tem- 
peratures all during the measurement stroke of the 
bell. An average of measurements taken at different 



times and even places in the bell may not be suffi- 
ciently accurate if bell and gas temperatures are not 
equal. 

The third difference between the modified and 
ordinary prover is in the location and design of the 
diverter valve. It is placed on top of the bell instead 
of on that portion of the inlet pipe that is external 
to the prover's tank. It is proposed that it be an 
electro-mechanical device designed to be able to 
remotely select and replicate the stroke time and 
flow rate history of the gas through the valve. A 
lateral discharge of the gas from four ports spaced 
90° around the new valve will be necessary to 
avoid a variable and disturbing force on the bell 
during its upward motion. A rotatable collar or 
other device with equally spaced and appropriately 
sized ports together with, possibly, a spring and/or 
an electric and/or a pneumatic motor to rotate the 
collar, will be needed to complete the essential 
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parts of the valve. It's not necessary that the new 
valve be located as shown in figure 1 (a), but at that 
location the initial balance and start-up procedures 
do not need to differ from that ordinarily used. If it 
is put on the external pipe of the modified prover 
some small changes in the initial balance procedure 
and in the gas flow equations would be needed. 

The differential equation of motion for the bell 
from reference [1], together with modified equa- 
tions of motion for the gas and sealing liquid, will 
be used to provide the basis to derive and compare 
the performance of the two provers. 

Special integrals of the Navier-Stokes (NS) equa- 
tion will be used (for both steady and unsteady mo- 
tion of the bell in the sealing liquid) to derive the 
liquid drag forces for use in the bell's equation of 
motion. Forces connected with entrained liquid 
and effects not included in the special integrals will 
be estimated for the same application. The NS 
equation also provides a justification to treat the 
motion of the sealing liquid as if it is a "solid" and 
thus a one-dimensional, second-order differential 
equation is used for the liquid motion to calculate 
its effect on the gas volume under the bell. 

Each set of differential equations will be sub- 
jected to numerical integration using the Runge- 
Kutta fourth-order method outlined in [2]. 



2. Equations of Motion 

2,1 Gas Motion 

The flow measurement process is initiated by a 
gradual closing of the diverter valve to start the 
upward motion of the bell from its balanced posi- 
tion. A calculation of the net mass flow into the 
bell during the stroke of the diverter valve is based 
on the difference between the initial steady state 
mass flow out of the valve and the smaller flow out 
during its stroke. If the ports are rectangular and if 
the stroke is assumed to be one of constant velocity 
the orifice port area A^ then can be calculated as 
A^=A^i (1— y). Here y is used as t/t^, the time / 
from start of the stroke divided by the total stroke 
time fv. and A^-, is the initial open area. For a con- 
stant acceleration stroke, A^=A^i (1— y^)- The lat- 
ter assumption is used in the subsequent 
calculations of the flow as its accomplishment 
seemed to present a simpler valve design problem, 
and possibly because it would provide a 
"smoother" acceleration of the bell upward from 
its float position. A valve might be designed for a 
still smoother operation to close slowly at both the 
beginning and end of its stroke, with a port open 



area represented, for example, by a function such as 
^v=^vi (l/2-(l/4) cosV'y + (3/4) cos ttj). The 
rate of area change with a valve so designed would 
be =0 at the beginning and end. Effects of valve 
design on prover performance probably will be re- 
searched better by calculations than by experiment. 

2.2 The Ordinary Prover 

The net mass flow rate into the bell can be writ- 
ten, using the well known equation for the orifice 
with coefficient of discharge Co, as 



dm/d/ = CD^vi(2piA/'i)'^' 
-Co^,i(2pA/')'^^(l-y^) 



(1) 



where AP is {P -P^ and APj is {P-P^ and the first 
term on the right side of eq (1) is the input gas 
flowrate FRI. This rate is to be compared to that 
established by the calculations from the derived 
motion of the bell and sealing liquid and for the 
change of state of the gas in the bell. The mass of 
gas in the bell and inlet pipe is represented by m 
and p is its density. 

The gas law is now used to relate P to the bell 
and oil motions so that it can be inserted into eq (1). 
The gas law can be written, using P^ as the initial 
balance pressure, as 



{P/P^ = {mV,/m,V)'' 



(2) 



where m, and F, are the initial mass and gas volume 
under the bell and inlet pipe, and Ax is the volume 
change under the bell due to its motion x from 
x—Q. AFis the change of gas volume from oil en- 
trainment on the bell as it moves upward plus the 
volume change caused by motion of the liquid seal. 
The exponent k may be taken as the ratio of gas 
specific heats when an apparatus is used that would 
give adiabatic conditions. For a very small appara- 
tus in which the process might be one of almost 
constant temperature, k should approach unity. 
Note here that A is the effective area of a bell 
prover and is given in [1] as 



^=^1+^2^3/(^3+^4) 



(3) 



where Ai is the interior cross section area of the 
bell, A2 the cross section area of the metal of the 
bell, A-i and Ai, are the inside and outside horizontal 
surface areas of the sealing liquid, respectively. 
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Equation (2) is based first on an assumed unifor- 
mity of conditions throughout both volumes Vi and 
Vtf which allows use of the heat-exchanger but no 
metallic-wool or porous-plug. Second, the fre- 
quency of a Helmholtz resonator built like the bell 
and inlet pipe was calculated to range from about 8 
to 3 times the expected frequency of the pressure 
pulsations as expected from the motions of the bell. 
Resonance at the Helmholtz frequency probably 
will not be observed experimentally. 

The gas volume correction for the entrained oil 
can be derived from eq (9) in reference [1]' assum- 
ing constant pressure (and density) during motion 
of the bell. Second, if the bell is held stationary 
while the pressure is changed, the oil seal will 
move and this also will change the volume of gas 
under the bell. It is not assumed here that the oil 
level change will necessarily be such as to maintain 
equilibrium with the pressure change. Thus A F in 
eq (2) can be written 

A V= [(A,-A,)/iA,+A,)] X Af-A,(h,-h,,{) (4) 

n 

where the summation indicates the progressive ef- 
fect starting from the initial balance condition for 
the bell. The entrained oil volume 2 A/ is the vol- 
ume on one side of the bell, and (hi—hii) represents 
the change of sealing liquid level from the initial 
balance position. 

Equation (19a) of [1] gives 2 A/ (assumed equal 
to 2A/3 and equal to DAy^) as 



2A/=(47r/(,)'^^ (Z3) Vo°'^ ;£ x'-^Af 



(5) 



where Zj is a dimensional constant and n is a coun- 
ter for the number of time intervals in the proposed 
summation process to replace the integral of eq 
(19a). Here x is the velocity of the bell, and Vq is the 
kinematic viscosity of the sealing oil. x has been 
used in the summation process assuming it 
adequately represents the quantity (x—hi) in eq 
(19a). 

The gas law is used to express the gas density p 
in terms of its pressure in the bell and inlet pipe and 
then the flow equation (1) can be written 



' Note that the last term of eq (9) in reference [1] has a 
typographical error. A ratio of densities, not a product, is cor- 
rect. 



dm 
dt' 



--FRI-CuA^, [lpA{P/Pd''' {P/P. 



pjp:)r'{\-f). 



(6) 



It should be noted that Co will be specified as equal 
to zero at the instant that 7= 1 to simulate a closed 
diverter valve. Equations (4) and (5) can now be 
used to insert A V into eq (2) for P/P\ which will be 
used in eq (6) above. These three equations contain 
numerous terms to be derived subsequently in the 
discussions of the liquid motions and the equations 
of motion for the bell. 

2.3 The Modified Prover 

A computer simulation of the modified prover's 
performance requires that eq (6) be replaced by 
two flow rate equations. The first is an equation for 
the net rate of accumulation of mass in V^ and is 
expressed as 



dm^/dt =FRI-FRI(P^-P)/(P^i-Pd 



(6a) 



where the second right-hand term represents both 
the mass rate of flow from the porous plug and also 
the mass rate entering the space under the bell. 
Second, the net rate of accumulation of mass m 
under the bell is 



dm /dt =FRI{P^-Py{P^, -P^ 
-CrA,i2pLPY\\-f) 



(6b) 



and again Co should be made equal to zero at the 
instant y is specified to equal unity. Note that 
(Pbi—P^ is the initial pressure drop that is to be 
specified across the porous plug and that the flow 
rate through it is proportional to (Pb—P). This re- 
quires a plug designed for viscous flow which 
probably would be more effective than one de- 
signed for turbulent flow. 

An eq (2a) can now be written similar to eq (2) to 
derive P as expressed by the equation 



P/Pi=(,m/mif (V/iVi+Ax + AV)'' 



(2a) 



with the reminder that P, m, m, are for conditions 
under the bell only. For the specified constant gas 
temperature in V^, P^ can be written as 



Pb=-Pbi'«b/'«bi- 



(2b) 



Equations (4) and (5) are used again in the equa- 
tions above as in those equations for the ordinary 
prover. 
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2.4 Liquid Motions and Forces 

It is intended here to derive the forces on the bell 
that are associated with the motion of the bell in 
the sealing liquid. For this purpose some aspects of 
Couette flow in both steady and unsteady state mo- 
tion between parallel plates will be discussed. It is 
also necessary to look at the development of the 
laminar boundary layer on a flat plate in the event 
that this could be dominant as to the force it applies 
to the bell. These considerations might be useful in 
examining further the forces from oscillatory liquid 
motion, but these will be ignored as they tend to 
cancel when both sides of the bell are considered. 
However, the liquid oscillation would affect the 
volume of gas under the bell and an attempt will be 
made to account for this effect. Other force aspects 
of the liquid motion, such as from entrained oil on 
the sides of the bell and motion of the oil at its 
bottom edge are estimated to include them in the 
bell's motion equation. 

2.4.1 Couette Flow I The well-known mo- 
mentum equation of NS may be used to derive the 
velocity profile for a fully-developed, steady flow 
between parallel plates. This condition is practi- 
cally satisfied if, in figures la and lb, r > >2Z> and 
(L— x)>>2Z>. These length conditions would 
limit the influence of the complex flows that exist 
at the free surface and at the bottom edge of the 
bell. Thus the applicable equation is written as 

— aP/poSx -I- vo (a^w/ax^-l- a^M/a;;^-l- a^M/az^) 



= M aw/ax -|--oaw/a;; -l-aw/a/ 



(7) 



where y is the distance from the bell's surface, z is 
along its circumference and u and v are the liquid 
velocities in the x- and ;> -directions, respectively, 
as illustrated in figure lb. For the specified condi- 
tions, the terms aw Vax^ d^u/dz^, du/dx, v and aw/ 
dt in eq (7) are all equal to zero. Two boundary 
conditions, namely that u = U at y=0 (or that the 
liquid on the bell's surface is carried up with the 
bell whose velocity is U) and second that w =0 at 
y—2b are now used in the integration of eq (7) 
which gives the profile as follows: 

u/U =\-{\/lbU){U + [lb^/lia]dP/dx)y 

+ {\/llioU){dP/dx)y\ (8) 

It should be noted that dP/dx was assumed to be 
constant across y in the derivation of eq (8). 



Equation (8) may now be used to specify a value 
of dP/dx such that the liquid up-flow adjacent to 
the bell is balanced by a down-flow adjacent to the 
tank wall. Thus, an integration of the flowrate 
across y derives the volume rate total which, when 
set to zero, gives dP/dx as 



dP/dx = 3iiU/2b\ 



(9) 



The pressure drop along the up-flow as expressed 
in terms of a liquid head Ah is derived from eq (9) 
to be equal to 'i\xU{L—xyib^p^. The up-flow 
must spill over at the free surface to build a liquid 
head at the tank wall comparable to Ah to provide 
for the flow down along the tank wall and for its 
recirculation at the bottom of the bell. When the 
value of zP/Zx is inserted into eq (8) the velocity 
profile for zero net flow is found to be 

w/i7=l+ (3/4^-27}, '(]=y/b and ^=0. (10) 

Some of these concepts concerning the motion of 
the sealing liquid are illustrated in figure lb. 

Several interesting features of the flow can now 
be derived from the velocity profile. First, the ve- 
locity gradient at the bell's surface is 



aw/a;; = -2i7/Z>,;;=0and 0=0 



(11) 



for calculation of the shear stress there. Second, 
the volume rate of flow (spillover) carried up by 
movement of the bell, between the surface and the 
w =0 point (v=2Z>/3), is calculated to be equal to 



= O.296/Z>i7,;; = Oto;;=2Z>/3 



(12) 



for one side of the bell. The quantity / is an aver- 
age circumference of the bell. Third, this volume 
rate can be compared to the entrainment rate d//d/ 
of liquid that adheres to the bell as its surface 
leaves the liquid. Reference [1] gives this rate for a 
constant bell velocity, x, as 



Af/At^Z^ev^'^x 



0.644 '1.644 



(13) 



For this prover, the entrainment rate ranges from 
1.4 to 6 percent of Q in the range of bell velocity 
researched and discussed in section 3. 

The equations of motion will be integrated nu- 
merically using a stepwise procedure between time 
intervals At. In a given interval the bell will have 
velocity i:= U which will be used to calculate force 
Fx on the bell during the interval. Equation (11) 
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gives the velocity gradient and when it is multi- 
plied by jLi and the wall areas involved gives F^ on 
the two sides of the bell as 



Fi=-4fx/(L-;c)C//6. 



(14) 



2.4.2 Couette Flow II (a) H. Schlichting [3] in 
his figure 5.6 and eq 5.24 presents the velocity pro- 
file in the fluid between a suddenly accelerated 
wall and one that is stationary. This "sudden" con- 
cept is regarded as applicable here to determine the 
force on the bell for the integration of the equation 
of its motion, because the bell velocity changes 
from interval to interval. His analysis started with 
two terms of the NS equation, namely the local 
acceleration term interacting with the friction 
force. This applies because the convective acceler- 
ation and pressure gradient terms are zero in this 
non-steady, parallel-flow situation. 

With that modified NS equation and using ap- 
propriate boundary conditions the derived velocity 
profile is given as a series of complementary error 
functions (erfc) of the quantities ai [=2i/2(vo?)'^^] 
and a [=y/2(voty^^]. Partial differentiation (with 
respect to j;) gives the velocity gradient at the 
moving wall O'=0) as 

du/dy = -(Ajc/(7rVoO'^')(l + 2e-<'°''' 

+ 2e-^'^'^\..),y=0 (15) 

where Ax represents the velocity "jump" caused 
by the "sudden" acceleration. The gradient— >co as 
t—*0, and as ?— >oo, the gradient decays to that for 
the steady moving wall. 

It is proposed here to calculate a time-averaged 
value of the velocity gradient from t=0 to t=At, 
and to use the "averaged" gradient to calculate 
force F2 on the bell. The time-averaged value cal- 
culates to be 

Cy^) = 2Ax/(7rvoAO '^' - 4Axe-*''^''''M73-VoAO'' 



r=0 to t = At 



(16) 



where the second and other terms on the right side 
of eq (16) are negligible for the small values of At 
to be used here. This result is taken as most 
applicable during the acceleration phase of the mo- 
tion of the bell. Thus eq (16) leads to 



i?2= -4ju-/^wHAA:/(7rvoAO'' 



(17) 



(b) The "sudden" concept is also regarded as 
applicable to estimate the motion of the oil seal (as 
caused by pressure changes in the bell) and the re- 
sultant change of h^ as used in eq (4) for the volume 
of gas under the bell. Here it is the oil that is moved 
suddenly, not the bell. If the velocity of the oil is 
constant over most of the space between the bell 
and the walls of the tank, a simple motion equation 
can be used to estimate A3 and hi by considering the 
oil as a solid in motion. 

A velocity profile in this "sudden" case is ade- 
quately represented (by a single term of a series 
again for the small time intervals to be used here) 
as 



f*2ai—a 
\l/2\ I --n2. 



u/hM^ATry'^) e-^^dTj 



(18) 



where the approximate equality sign also signifies 
that A3 is not exactly constant across the distance 
2b . The integration limits for tj are calculated from 
tti and a as given previously. A numerical integra- 
tion of eq (18) for values of t assumed equal to 
integration time intervals At to be used here does 
show an almost constant velocity across practically 
all of the distance 2b . Thus, a value of A3 derived 
from the "solid" equation of motion should be re- 
duced by a factor derived from an average value of 
u /A3 across the distance 2b . It should be noted here 
that the oil motion effect on the gas volume is ex- 
pected to have most significance while the diverter 
valve is closing. If the subsequent motion of the 
bell is at constant velocity the motion of the oil seal 
should have insignificant influence. Note the big 
"if in the previous sentence. 

2.4.3 Laminar Boundary Layer A. H. Shapiro 
[4] presents information on the development of this 
layer and its velocity profile which was used to 
estimate its thickness and possible drag force on the 
bell. It can be shown that the boundary layer thick- 
ness quickly becomes comparable to the separation 
between bell and tank. This would require large 
and unrealistic pressure gradients in the region 
such as to void the computations. These and other 
considerations lead to the conclusion that the lami- 
nar layer cannot build up in the apparatus discussed 
here. 

2.4.4 Other Forces (a) Sealing liquid that is 
entrained on the walls of the bell previously was 
stated to modify the indicated volume of gas col- 
lected. Here the liquid acts to modify the mass of 
the bell by way of the force of gravity on the liq- 
uid. This will be the third non-conservative force 
F3 calculated below: 
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= -2pog(47rA0"' (Z3) vo"*^ X^'-^-^A^ (19) 

n 

Ob) Liquid at the bottom edge of the bell is 
sucked upward by it which creates a drag force 
there which can be stated as 



F4= -A2poX^/2. 



2.5 Bell and Liquid Motions 



(20) 



The equation of motion (13) of [1] together with 
the initial conditions and static balance eqs (14), 
(15), and (16) of [1] will be used here to set up the 
second order motion equation for the bell. As indi- 
cated in section 2.4.2, it appears that the oil motion 
can be regarded as the motion of a solid and this 
basis will be used to set up a second-order differen- 
tial equation of motion for the oil. These equations, 
in conjunction with the first-order differential eq 
(6) for the gas motion, will comprise the system of 
equations to be handled by numerical integration. 
For this purpose they will be transformed by the 
usual procedure into five first-order equations, 

For the motion of the bell eq (14) of [1] was used 
first, here with x,=0, to calculate h^i at the initial 
balance condition using an assumed (or specified) 
value of (Pi-PJ. Second, eqs (15) and (16) of [1] 
were used to calculate Mj/Lj (the chain mass per 
unit of length) so that the cam and buoyancy com- 
pensator M4 could be replaced with the chain com- 
pensator. By this means the term involving x^ is 
removed from the bell's equation of motion. The 
equation thus will be a little simpler. Performance 
with the chain can be compared later with that of 
the usual suspension system. Third, eq (13) of [1] 
then can be simplified further to read as below (af- 
ter noting that Ax (Pi—P^)=g (Mx—Mi—poAiL) 
-pogAih^^: 



(Mx+M2+I/ri)x= -{p^2-2M^/L^)gx 

+^,(P-Pi) + pQg^2(/24-h4i) + 2F. 



(21) 



All terms above are zero at the initial balance con- 
ditions of x='LF=0 and h4=h4_\ and P=Pi. EF 
designates the sum of the drag forces. Figure 1 lo- 
cates and defines some of the quantities used. / is 
the moment of inertia of the pulley. 

Two different approaches might be selected to 
set up the oil motion equation, one regarding the 
oils on the two sides of the bell as two separate 



"solids", or by treating the oil as a single "solid" 
using the assumption that ^3=^4. For the two 
"solids" case the equation is quite complicated and 
contains a product of the oil acceleration and ele- 
vation and also its velocity squared. This seems an 
unwarranted complication at the present time and 
the single "solid" will be used. 

The equation of motion for the single "solid" is 
derived using the familiar concept that acceleration 
is given by force per unit mass and that 
(^4-^41)=- (^3-^31), to give 

'li,=(P-Pd/[2po(L -x)]-g(h,-h,dAL -X) 

-Force Loss/|>o(-43-l-A)(^-^)]- (22) 

It should be noted here that the oil mass to be 
moved is poiA^+A^XL —x), and that the term con- 
taining /J4 is the liquid's back pressure. The factor 2 
(and also unity in the liquid head term) arises 
from the pressures acting on the assumed 
area=(l/2)(^3-|-^4). A "Force Loss" is included in 
the equation to include the effect of the turn that 
the oil must make at the bottom of the bell of the 
ordinary prover. This loss is assumed to be equal to 
one-half of the velocity pressure, or (l/2)(l/2) pohj 
multiplied by the area A^, which then gives the last 
term on the right side of eq (22) above as (1/8)^4 
I A4 1 /(L —x). The quantity | A4 1 is used to make the 
force loss applicable to both flow directions. When 
the two orifices are installed to constrict the oil 
flow space down to an area of (1/2)^5 on each side 
of the bell to modify this part of the prover, a pres- 
sure drop will be introduced into the flow there. 
This drop will arise from the frictional loss of me- 
chanical energy mostly on the downstream side of 
each constriction. For an area constriction down to 
about 10 percent, a possibly conservative estimate 
of the pressure drop would correspond to one ve- 
locity head pressure based on the liquid's average 
velocity in the constricted area. On this bases the 
last term of eq (22) will be equal to the quantity 
(l/2)(^3-|-^4)^ h^lh^l/AjiL—x) when the orifices 
are used. 

In the context of the discussion above it should 
be noted that the "solid" oil motion will have some 
effect on the ordinary bell by the necessity of mak- 
ing the U-turn there. If it's assumed that the down- 
ward motion of the oil pressurizes the "dead" oil 
space below the bell and acts to help move the bell, 
the force Fs can be estimated as 



F, = (m)poA2hl 



(23) 
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where the above is based on one-half of the veloc- 
ity pressure. The above force is replaced, for the 
modified bell, by a suction force on the bottom of 
the bell as caused by the downward motion of the 
oil through each orifice constriction. This motion is 
induced by the upward motion of the bell. As dis- 
cussed above for a proposed area constriction in 
the orifice down to 10 percent, one velocity head 
pressure loss (based on the velocity in the constric- 
tion) is probably a conservative estimate. On this 
basis Fs for the modified prover is 



F,= -(\n)p^2iA2/A,fx\ 



(23a) 



These forces F5 are introduced as reminders that 
they and other "end effect" forces exist but proba- 
bly have small effects. 

2.6 Summary of Motion Equations 

The gas flow equation (6) for the ordinary 
prover is 

m=FRI-CoA,, (IpAY^ (P/PyY'' (P/P, 

-PJP,y^'(l-y') (24) 

where Cd = for the condition y>l and where the 
factor (1— 7^) assumes the use of a constant accel- 
eration diverter valve. Equation (24) needs auxil- 
iary eq (2) for P/Pi and eqs (4) and (5) for V/V^ as 
follows: 



P/Pi = [m/mt/{V/Vi)Y; 



V/V,= l+Ax/V, + (Z,/V02(x''^^t) 



(25) 



+Z2(A,/Vdh,; (26) 

where Z, = [(A,-A,)/{A,+A,)\ (Z3) (477^,)"' ^o"'^; 

Z2 = a correction factor for lack of oil "solidity"; 
Zi = a constant whose value depends on the sys- 
tem of units that will be used. 

Note above that A4 has replaced (ha — hi) because 
Aj is assumed equal, in eq (26), to A4 and that the 
motion equations when integrated will give 
h4=(h4-h^;). 

The gas flow eqs (6a) and (6b) for the modified 
prover are: 

m,=FRI-FRI {P,/p,-P/p;)/{P,yp.,-\); (24a) 



and 



m =FRI(,P^P-P/P^/(,P^,/P- 1) 

-C^A2p,PiP/P,)'\P/P-PJP]''\\-y\ 

(24b) 

Equations (24a) and (24b) need auxiliary eqs (2a) 
and (2b) for P/P, and /"b/P,-, respectively, and eqs 
(4) and (5) also will be used again in the following 
equations: 



P^P; = {Pu/P^m^m^, 
P/P; = [m/my(V/V,)f 



(25a) 
(25b) 



V/K=l+Ax/V;+(Z,/V{)'2(x'''*At) 

+ Z2(,A,/Vdh, (26a) 

The equation for the motion of the bell eq (21) is 

(M, + M2 + lM)x=A,(P-Pi) 

5 
-{poA2-2M,/L,)gx+pogA2h,+J^F (27) 

1 

which uses individual drag forces, as follows: 
F, = - (4jLt f/b)(L -x)x; (27a) 

F2=- (4j[A ^/(TTVoy^XL -x)Lx/{Lty'^- (27b) 

^3= -2pa?(47r^,)"' (Z3) v^'"'J^{x''^Lty, (27c) 

n 

F4=-(l/2)poA2x'; (27d) 

F5 = (1/2) poAzhl (ordinary prover) (27e) 



F,= -{\/2)poA2{,A2/A,yx' 
prover). 



(modified 
(27f) 



The equation of motion eq (22) for the sealing oil in 
the ordinary prover is 

h=i,P-P^/[lpo(L -x)]-gh4/{L -X) 

-(\/8)k\k\/(L-x) (28) 

and for the modified prover the last term on the 
right side of eq (28) is specified as equal to 
50A4|A4(L —x) when A^ is constructed to be equal 
to 0.1 (A}+A4). With this substitute, last term of 
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eq (28) in place we may label the motion equation 
for the sealing oil as eq (28a). 

Equations (24), (27), and (28) or their (a) and (b) 
substitutes for the modified prover, comprise the 
two systems of equations to be integrated numeri- 
cally by the fourth-order Runge-Kutta method as 
outlined in [2]. 

A recapitulation of the main assumptions used in 
the derivations of the above equations seems advis- 
able. It does not seem advisable at the present time 
to develop a change to eq (24) to account for a 
possible non-constant value of Co. Equations (27b) 
and (28), and the far right-hand term of eq (26) are 
based on a special integral of the NS equation for 
the case of "sudden" motion of the bell and of the 
sealing liquid. Equation (27a) was derived from the 
NS eq (7) for steady flow between parallel plates 
without consideration of end effects. Equation 
(27e) and the far right-hand term of eq (28) but not 
for eq (28a) are used as reminders of end effects not 
fully accounted for. The forces of these "esti- 
mated" end effects turned out to have an insignifi- 
cant effect on the calculated performance of the 



prover. Except in eqs (27), (26a) and (3), A^ was 
assumed equal to A4 with corresponding assump- 
tions /3=/4 as in eq (27c) and —hi— ha, as in eq (28). 
And, finally, the flow rate into the bell of the ordi- 
nary prover is assumed to be constant and unaf- 
fected by the probable small changes of state of the 
gas in Fb. 



3. Calculated Results 

The numerical integrations of the equations of 
motion are used to derive the performance of the 
provers which will be illustrated by various curves 
showing the motions of the bell and sealing oil and 
the derived gas pressure and gas volume correction 
all expressed as a function of the time. (See figs 2 to 
5). These figures can give a good qualitative illus- 
tration of the performance, but further computa- 
tions are necessary to derive quantitative 
evaluations and comparisons of the measurement 
performance. Equation (2) can be rearranged for 
this purpose so that the true mass collected can be 




.99900 



01234567 

Time(s) 

Figure 2. Bell pressure, and height of sealing liquid, plotted vs time. Program HIFLO. 
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Figure 3. Actual volume of gas collected, and volume correction, plotted vs time. Program HIFLO. 



expressed in terms oi p-A (;C2— ;ci) with corrections 
applied to it for the fluctuations of P and A V. Sub- 
scripts 1 and 2 are used on ;c to signify values at the 
beginning and end of the chosen timing interval, 
respectively. It seemed preferable to separate the P 
and A V corrections, and to do this the pressure ra- 
tio in eq (2) was expressed as very closely equal to 
(l-f-(i'-i'i)Ai'i). With some use of algebra, eq (2) 
can be written for the modified prover in the form 



{m2-m^=^pA{x2-x^ 



1-H 



^(^2-^1) 



+ 



{^P^-P^{y,^Ax^ _ (,P,-PW:+Ax,^ 



■V 



kP-Aixi—x^ 
(Pki-Pu)Vb 



kPy4 (X2 — X1) 



{k = \)P,^(X2-X) 



(29) 



and when the ordinary prover is used Fb = in the 
far right-hand term and F, should be replaced by F, 
which is given by (F-f- V^). 



3.1 The Ordinary Prover 

Almost complete physical data were available to 
the authors on a bell prover of nominal capacity Vi 
of 5 ft' (0.1415851 m') and these were used in three 
computer runs for the first look at the utility of the 
proposed analytical method. Results of these runs 
indicate the effectiveness of the integration method 
as well as demonstrate the size and sources of er- 
rors likely to be encountered with the usual labora- 
tory procedures. The prover specifications are 
listed in table 1 together with data on initial bal- 
ance conditions, fluid properties and the computer 
run conditions used in the programs HIFLO and 
LOFLO. Nominal measured stroke times are 5 and 
50 s, respectively, for these two programs. 

The two runs at the high rate of flow (nominal 
collection time of 5 s) were done with different 
time intervals in the integration process, namely 
0.002 and 0.01 s. Rates of flow derived from use of 
these two integration processes differed by about 
0.01 percent. Comparison of other physical vari- 
ables also showed differences near 0.01 percent. On 
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the bases of this good agreement it was decided to 
use the 0.01-s time interval for the low (50 s collec- 
tion) flow rate calculation. It seemed reasonable to 
expect much larger differences than 0.01 percent 
between prover performance at the two rates of 
flow. 

Figures (2) and (3) are graphical representations 
of the calculated results at the high flow rate and 
figure (4) presents results at the low flow rate. Fig- 
ure (2) shows a relatively large initial fluctuation of 
the pressure, about ±46 percent of APi from P^ 



with a relatively long period, initially about 1.8 s 
and decreasing to about 1.1s. The pressure fluctua- 
tion decreases to about ±30 percent of APi from 
the mean pressure, but the curve near 7 s shows a 
hint of a possible increase at that point. Higher fre- 
quency fluctuations of the pressure, with periods of 
about 0.3 s, are also evident. The amplitude of the 
oil motion as converted from values of hi to pres- 
sure is about 120 percent of the pressure fluctua- 
tion. It is evident that the pressure fluctuations are 
strongly modulated by the oil fluctuations. It is also 



Table 1. Prover specifications and computer run conditions 



Prover specifications' 


A 


A, 




Ai 




A, 




A, 


/ 


.2238913 
2.409936 


.2183231 
2.35 




.01189164 
.128 


.05481303 
.59 




.06224530 1.678537 SI 
.67 5.507 U.S. 


M' 


M, 




M, 




M^/Lf 


M, 




I/rt 


94.71400 
6.49 


99.23840 
6.8 


85.45006 
5.855197 




5.569795 
.1163280 








5.032009 SI 
.3443026 U.S. 


b 


Kb' 


^5/(^3+^4) 


K, 










0.0175 
0.0575 


0.0062 
0.22 


0.1 
0.1 




0.14 m 
5 ft' 


' SI 
U.S. 








Initial conditions and fluid properties 


L" 


g 




P. 




(^i-^a) 




A4i 


X-' 


.82296170 
2.7 


9.8054 
32.17 




101362.1 
2117 




239.4003 
5 


1.344412(10 
4.410790(10- 


2) SI 
-^) U.S. 


V^ 




OTi 




k' 


Vo 






Ho 


2.208728(10-^) 
0.78 


3.409868(10-^) 
2.336505(10-') 




1.4 
1.4 


.2006714(10 
2.16(10-4) 


-'■) 




170.453(10-") SI 
3.56(10-*) U.S. 


Po 


Pa 




Pi 




(Pb,-PO' 








850.3685 
1.65 


1.201339 
2.331(10-') 




1.204176 
2.336505(10-') 


40(Pi-P,) 

40(Pi-Pa) 


SI 
U.S. 






Computer run conditions 


Program name 


m=FRI 




fv 


Z2 


M 




Nominal Collection 
60OT/pa time, pa V\/m 



'\ 



HIFLO" 

HOLHOG 

LOFLO 



3.4018341(10-^) 


2 


.95 


.002 


2.331(10-') 


2 


.95 


.002 


Same as above 


.2 


.95 


.0002 


3.4018341(10-') 


2 


.90 


.01 


2.331(10") 


2 


.90 


.01 



1.699021 5 s SI 

60 scfm 5 s U.S. 

— — same as above— 

.1699021 50 s SI 

6 scfm 50 s U.S. 



' SI units are kilogram, meter, newton and seconds (MKS). U.S. units are slug, foot, pound-force, seconds. 
' Program HIFLO was also done with an integration time interval At of 0.01 s, 
= Derived or assumed as probably applicable. 
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Figure 4. Bell pressure, and gas volume correction, plotted vs time. Program LOFLO. 



strongly evident from figures 2 and 3 that the gas 
volume correction AV under the bell fluctuates 
practically in unison with h4. This fact suggests that 
failure to measure h^ possibly can be the major 
source of error. 

At the low rate of flow figure 4 shows a much 
smaller initial fluctuation of pressure, about ±4.0 
percent of APj from Pj, and again a relatively long 
period initially of about 1.8 s which decreases to 
about 1.4 s. The pressure fluctuation and the pres- 
sure level both decrease, the former to ±1.3 per- 
cent of APi about the mean pressure which is about 
9 percent of AP, below Pj. These fluctuation ampli- 
tudes are about six percent of those at the high rate 
of flow, as compared to the flowrate ratio of 10 
percent. The amplitude of the sealing liquid motion 
(not shown) begins at about 115 percent of the 
pressure, but winds up near the end of the measure- 
ment stroke at about twice that of the pressure fluc- 
tuation. As at the higher flow rate, the gas volume 
correction (not ordinarily measured) fluctuates 



practically in unison with h4 and with the low fre- 
quency pressure fluctuations. The high-frequency 
fluctuations have periods of about 0. 1 9 s which in- 
creases to about 0.31 s. 

The high frequency motion calculated for the 
bell can be shown to be real by considering the bell 
to be vibrating as if suspended from a "spring" 
which can "replace" the air enclosed in the bell. 
For the case of the very low values of the pressure 
fluctuations prevailing here, this concept is suffi- 
ciently accurate to derive the "spring constant" on 
X as equal to kAiPyV where Kis the volume of air 
and X is the motion of the bell about a stationary 
position. The vibration period for this well known 
case is iTr/ikAfP/iViM.+Mj+lM)]}'^^ which 
calculates to be in the range from 0.2 to 0.4 s. Inclu- 
sion of the buoyancy force and damping (drag) 
forces would not change the calculated period suf- 
ficiently to warrant the added complication. Here 
again these periods are in conformity with those 
derived from the full motion analysis. 
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Figure 5. Bell pressure, and gas volume correction, plotted vs time. Program HOLHOG. 



An inspection of eq (29) and figures 2, 3, and 4 
demonstrates the problems involved with making 
accurate measurements of the fluctuating physical 
variables and also the need to correlate those mea- 
surements with the timing points used to define the 
stroke of the bell. It is also evident that various 
choices and measurement methods can be used to 
try to derive the various corrections in eq (29) for 
the true collected mass. 

3.1.1 Method 1 If the measurements of P and 
either Lhi, (or A V) are accurate at any instant and 
accurately correlated with the bell positions X] and 
Xi, eq (29) provides the true mass. 

3.1.2 Method 2 Figure 2 demonstrates that no 
corrections to p^ixi—x^ are necessary if the ini- 
tial and final timing points (at 3.36 and 6.2 s, for 
HIFLO) can be identified by the instrumentation 
when the variables of (P—Pi) and Ah^ both equal to 
zero at both bell positions (a zero-zero criterion). 
In the more general case (LOFLO, for example) 
this criterion does not prevail. However, pressure 
transducers probably are available with capability 
to follow the low frequency fluctuations of the 



pressure (as caused by the motion of the oil) with 
good accuracy especially at the fluctuation mid- 
points. If the timing points Xj and X2 are made to 
coincide in time with those mid-point pressures, the 
corrections derived with eq (29) would also be ac- 
curate because the mid-points of the A Fand the AP 
motions are shown here to coincide. These meth- 
ods should be used but apparently are not because 
of the practical difficulties. A third procedural 
compromise. Method 3, is usually used to derive 
the corrections. 

3.1.3 Method 3 Measurements are made with 
a water manometer to give (what is believed to be) 
mid-point values of the high-frequency fluctuating 
pressures (AFP) at two preselected, convenient po- 
sitions Xi and Xi- These measured AFP values are 
used in eq (29) along with calculated values of A F 
(based on AFP) to derive "corrections" to 
p^(x2—Xi). The "worst" timing point choices that 
can be derived from figures 2 and 4 are used to 
illustrate in table 2 the errors that might result with 
use of Method 3. "True" corrections (of Method 1) 
are also listed for comparison in Table 2. 
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Table 2. Measurement corrections for worst timing choices 



Program 



Correction 
Method 



A(AK)% A(AP)% 



APb% 



% Total, 
correction 



HIFLO 1, True -0.752 -0.072 

HIFLO 3, AFP -0.196 -0.065 

HIFLO Remaining error after AFP Correction = 

LOFLO 1, True -0.017 -0.060 

LOFLO 3, AFP -0.053 -0.001 

LOFLO Remaining error after AFP Correction = 

HOLHOG 1, True -0.138 -0.045 

HOLHOG 3, AFP -0.108 -0.022 

HOLHOG Remaining error after AFP Correction = 



NA 


-0.824 


3.73 


6.95 


NA 


-0.261 


3.73 


6.95 


0.563 








NA 


-0.077 


3.8 


44.6 


NA 


-0.054 


3.8 


44.6 


0.023% 








0.001 


-0.184 


3.73 


6.6 


0.001 


-0.131 


3.73 


6.6 


0.053% 









For program HIFLO, one "worst" timing point 
choice would need a very large true correction, 
while the correction based on values of the AFP is 
only 32 percent of what is actually needed. For a 
LOFLO "worst" choice, use of the AFP gives a 
correction that is 70 percent of the true correction. 
The net result is that flow rates based on AFP val- 
ues would be higher than actual by 0.56 and 0.02 
percent for the flow rates HIFLO and LOFLO, 
respectively. It should be noted the corrections for 
entrained oil on the sides of the bell amounted to 
—0.013 and —0.003 percent for those two rates. 
These are not shown above because they are in- 
cluded in the values of A F in the true correction 
procedure [see eq (4)]. It should also be noted that 
other "worst choices" would give corrections of 
opposite sign. 

There are three suggestions to be made concern- 
ing possibilities for improving the instrumentation 
and/or measurement procedures used with the or- 
dinary prover. 

3.1.4 Suggestion 1 Instrumentation problems 
would probably prevent implementation of the 
ideal situation in which Method 1 could be effec- 
tively used. However, experiments should be con- 
ducted to try to overcome these difficulties. 

3.1.5 Suggestion 2 It might be possible to im- 
plement the zero-zero criterion as described in 
method 2 to use what is essentially a "correction- 
less" procedure. Two requirements are necessary 
for this purpose. First, the prover must be modified 
to maintain the average of the maximum/minimum 
pressures, P^, constant and equal to Pi during the 
stroke of the bell. Second, an instrument must be 
available to be able to identify bell positions where 
the fluctuating values of AP and A/z4 are equal to 
zero. The first requirement probably can be met by 
an adjustment of the buoyancy weights (downward 
by about one and seven percent, respectively, as 
computed from the pressure curves for HIFLO 



and LOFLO). This should be and is confirmed if 
some bell positions, particularly near X2, can be 
found where AFP values or transducer measured 
values are equal to P^. These positions need not 
necessarily be those that are finally used for the 
rate measurements. 

An instrument which might be used for phase 
indications of AP and Ah4 is the optical encoder 
used by G. Kulin, P. Huang, and G. Mattingly [5] 
to better measure the volumetric displacement of 
the bell and its velocity. It was mounted on the 
pulley and it gave (almost) linear vertical motion in 
the form of electrical pulses. Their temporal traces 
of the encoder's output voltage displayed fluctua- 
tions of the bell's velocity at approximately the 
same period of about 1.8 s as calculated herein. 
Some computations for program LOFLO (not 
shown) demonstrated that maximum bell velocity 
occurred when both AFand values of (P— p) both 
were zero. Therefore, presumably, the encoder's 
indication of the maximum bell velocity also occurs 
at the zero-zero condition. This leads to the sugges- 
tion that the timing points should be made to coin- 
cide with the maximum velocity indications of the 
optical encoder to try to eliminate the otherwise 
needed, but in accurate, AFP corrections to derive 
the true mass of gas collected. As an alternative to 
the above, more attention should be given to fully 
implement Method 2. 

3.1.6 Suggestion 3 It appears to be possible to 
implement the zero-zero criterion as described in 
Method 2 using the AFP procedure of Method 3 
which can be modified to measure the true col- 
lected mass of gas. Two modifications appear of 
possible utility, both needing a constant value of 
Pn,=Pi during the stroke of the bell. The first 
would use computed motions of Ah4 (as described 
herein) to identify and use the bell positions Xi and 
X2 where AA4 values are expected to be zero. AFP 
values would then also equal P and corrections to 
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PiA (X2—X1) would also equal zero. Results from use 
of this method could then be compared to results 
obtained below with the second modified AFP pro- 
cedure. 

The second modified AFP procedure would uti- 
lize a multiplicity of strokes of the bell each with 
positions Xi and X2 systematically moved until both 
give AFP values equal to Pi. P^ must also be main- 
tained equal to P, with this method (see Suggestion 
2). This method is based on the assumption that 
AFP values of pressure fluctuate about Pj with a 
smaller amplitude than those of the high frequency 
pressure fluctuations, and with a period of the low 
frequency fluctuations of A/14. It's expected that on 
average five different strokes of the bell would be 
needed with Xi and X2 each moved each time a dis- 
tance corresponding to the distance traveled by the 
bell during 0.05 period of the Ah4 fluctuation. It's 
assumed in the experimental situation that the pres- 
sure measuring instrument actually does indicate 
the AFP values. If all of these criteria prevail the 
modified, multiple AFP method is also correction 
less. An advantage of these two procedure modifi- 
cations is that neither requires apparatus changes. 

3.1.7 The Modified Prover Figure 5 presents a 
graphical picture of the calculated temporal traces 
for the pressure and gas volume correction as 
derived in program HOLHOG which was speci- 
fied to have the same flow rate as in HIFLO. These 
results show that the modifications made to the 
prover caused large reductions in the amplitude of 
the pressure and oU motions (particularly for the 
latter) as was intended. Each oil orifice had a 10 
percent open area, the drop across the Porous Plug 
was specified as equal to 40A/'j, and t^ was reduced 
to 0.2 s (from 2 s as used in HIFLO and LOFLO). 
These three changes were made simultaneously 
and "almost" arbitrarily, and hence the name 
HOLHOG was chosen. Note that A V increased by 
only about 50 percent over that in HIFLO in spite 
of the 14-fold increase of AP at the end of 0.2 s. 
These larger upsurges, however are subsequently 
efficiently dampened. Although the frequencies of 
the motions of AF and A/* are the same as in the 
ordinary prover, the practically constant phase dif- 
ference between the two motions is not evident 
here. 

The calculated "worst" corrections for HOL- 
HOG listed in the table above demonstrate that the 
correction remaining and needed after application 
of the AFP Method 3 is gratifyingly reduced to 
0.053 percent from the 0.586 percent needed for 
the ordinary prover. This possible "worst" result 
for HOLHOG is improved to a needed remaining 
correction of only 0.006 percent if a "best" timing 



choice happens to be made. Phase differences be- 
tween the AFand AP motions that change during 
the stroke of the bell limit suggestions for improve- 
ment in methodology to Suggestion 3 above. How- 
ever, the need for improvement is much, much less 
than for HIFLO at its worst. 

It's appropriate at this point to discuss the "al- 
most" arbitrary group of apparatus changes (speci- 
fied for HOLHOG) and the problem of setting the 
initial conditions. As a preliminary to HOLHOG, 
computations of the effect of the orifice installation 
alone, and next of the effect of the porous plug 
alone, indicated that dampening was not effective 
enough before the bell reached X\ to make a signifi- 
cant improvement. A third preliminary calculation 
was done to evaluate the effect of a single change 
in an initial condition, namely a reduction of ^v 
from 2 to 0.2 s in program HIFLO. This change 
produced a two-fold increase in the initial upsurge 
of AF over the already large increase evident in 
HOLHOG. This large increase illustrates the im- 
portance of the initial condition, and second that 
the orifice effectively dampened the motion of the 
oil in HOLHOG and especially the initial upsurge 
motion. 

Other changes of initial conditions, such as in the 
flow rate characteristics of the diverter valve, and 
changes in the relationship of V-, to Vf, to F,, also 
could have significant effects on the performance 
of both provers studied here. This discussion 
should now revert back to the need to use a motor- 
ized and stroke-designed diverter valve as previ- 
ously suggested herein. This discussion also leads 
to the conclusion that prover performance cannot 
be adequately specified until the initial conditions 
are too. 



4. Conclusion 

It appears that the equations of motion for the 
moving elements of the prover system and the pro- 
cedure used in their numerical integration can pro- 
duce calculated results which offer a reasonable 
chance that they can be used to design improve- 
ments in the apparatus and measurement proce- 
dures. "Dynamic errors" in the measurement of 
flow rates can arise primarily from the fluctuating 
motions of the sealing liquid, second from fluctuat- 
ing velocity of the bell, and third and usually least 
in magnitude from entrainment of liquid on the 
bell. Proposed additions to the prover can be used 
to (1) assure needed equality between gas and bell 
temperatures and (2) to reduce the amplitude of the 
fluctuating motions as one way to reduce the mag- 
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nitude of measurement errors. Three other sugges- 
tions for improved measurement procedures in- 
clude an ideal situation in which pressures and oil 
motions are accurately measured at any two bell 
positions of choice. A second more practical sug- 
gestion would require that the timing points for the 
bell motion be properly correlated with indications 
of the mid-point pressures from a pressure trans- 
ducer or correlated with indications of the maxi- 
mums in the fluctuating bell velocity. A third 
suggestion for improved accuracy entails the uti- 
lization of two different approaches both being 
modifications of that currently in use which relies 
on manometer reading of the average indications of 
the high-frequency fluctuating pressures at the two 
bell positions of choice. Some changes in the initial 
conditions produced performance changes which 
emphasizes that blanket performance and accuracy 
statements cannot be specified unless all initial con- 
ditions are also specified. 



5. Nomenclature 



A-- 

Ar- 
A,- 

A,-. 
A,-- 
A,- 

A.-- 
AFP-- 

2b-- 
Cd= 

A = 

F-- 

FRI- 

/-- 

g-- 

M-- 
I-- 

k-- 



L: 

e-. 



-- See eq (3), effective bell area 

= Cross-sectional, open area of the bell 

= Horizontal cross-sectional area of the metal 

wall of the bell 

= Horizontal area of the inside sealing liquid 
= Horizontal area of the outside sealing liquid 
= Sum of orifice open areas in sealing liquid 

spaces 

= Area of ports of the diverter valve 
= An average or mid-point values of the fluc- 
tuations of P 

= Separation distance between bell and tank 
= Coefficient of discharge for the diverter 

ports 

= Change of 
= Force 

=Gas flowrate into the bell = mi 
= Volume of entrained oil on one side of the 

bell 

= Local acceleration of gravity 
= Sealing liquid elevations with reference to 

L 

=Oil "spillover head" 
= Moment of inertia of the pulley 
= An average ratio of isothermal to adiabatic 

compressibility — for air between 1.0 and 

1.4 
= Sealing liquid depth when a; =0 
= Length of the chain 
= Circumference of bell 



M= Mass of the bell 

Ml = Mass of the bell plus that portion of the 
chain which is on the same side of the pul- 
ley when a; = 

Mass of the counterweight plus its chain 
weight 

Mass of the chain 

Mass for buoyancy compensation (not used 
here) 

Mass of the gas in the bell 
Average of the maximum/minimum pres- 
sures in the bell 
Pressure in the bell 

Pressure in the bell inlet pipe of the modi- 
fied prover 
Differential pressure 

Volume rate of spillover liquid on one side 
of the bell 
seconds 
Time, seconds 

Diverter valve closing interval, seconds 
Integration time interval 
Time at beginning of timing interval 
Time at end of timing interval 
Bell velocity along x, or "freestream" ve- 
locity down along X 
Velocity of sealing liquid to positive x 
Initial volume of gas under the bell only 
Volume of gas in the bell including V\ or F, 
Volume of gas in the bell inlet pipe 
Nominal prover capacity 

Gas volume correction 
Velocity in the ^-direction 
Bell position at beginning of timing interval 
Bell position at end of timing interval 
Bell's movement upward from zero posi- 
tion 

Distance from bell's surface 
Distance around circumference of the bell 
A group of constants in eq (26a) 
A correction factor for lack of oil "solid- 

ity" 

1.09 m-°" s'-^' (SI) or 0.36 ft"" s'-^' (U.S.) 

j/[2(voO"'] 

2^/[2(voO"'] 

t/U 

y/b. Also, value of variable in eq (18) 

absolute viscosity of the sealing liquid, 

force second/length^ 

kinematic viscosity of the sealing liquid, 

lengthVsecond 

Liquid or gas density 



M2 = 

m = 

P= 

P,= 

AP= 
Q-- 

s = 
?= 

Ar= 
U = 

C/= 

M = 

Vr- 
V= 
V^- 
Vr- 
Vr- 
AF= 

D = 

^1 = 

x= 

y= 

z= 
Zr- 
Z,-- 

z,- 

a= 

ai = 

7= 

Vo= 
P = 
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Subscript 

a = Ambient condition or air 
b = conditions in inlet pipe 
i = Initial value, balanced condition 
n= Integration time interval counter 
o = Liquid 
v = Diverter valve 

Superscript 

' = First derivative with respect to time 
"= Second derivative with respect to time 

Dimensions Used 

slug (=14.593883 kg) 
foot (=0.3048009 m) 
Ibf (=4.4482290 N) 
Ibf/ft^ (=47.880056 N/m^) 
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